Neel order in doped quasi one-dimensional antiferromagnets 
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We study the Neel temperature of quasi one-dimensional S=l/2 antiferromagnets containing non- 
magnetic impurities. We first consider the temperature dependence of the staggered susceptibility 
of finite chains with open boundary conditions, which shows an interesting difference for even and 
odd length chains. We then use a mean field theory treatment to incorporate the three dimensional 
inter-chain couplings. The resulting Neel temperature shows a pronounced drop as a function of 
doping by up to a factor of 5. 
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The study of doped low-dimensional antiferromagnets 
has been a very active field since the discovery of high- 
T c superconductivity. A particularly simple form of 
doping results from replacing some magnetic Cu ions 
by non-magnetic ions like Zn. In this case the system 
is well described by the Heisenberg model with some 
spins removed from a regular lattice. This problem has 
been quite extensively studied both theoretically and ex- 
perimentally in the quasi-two-dimensional case. Recent 
Monte Carlo simulations have shown convincingly that 
the zero temperature 2 dimensional system remains Neel 
ordered for impurity concentrations, p up to the classical 
percolation threshold, (p c « .407, on the square lattice) 
jl], ^ . An analytic approach has been developed based on 
spin-wave theory and the T-matrix approximation, valid 
at impurity concentrations well below percolation, and 
extended to include weak inter-plane couplings |||. Re- 
cent experiments on La2Cui_ a; (Zn,Mg) 2 ;04 Q are consis- 
tent with the T = critical concentration corresponding 
to classical percolation, and agree in detail, at lower im- 
purity concentrations, with the analytic theory. 

Here we study the effect of doping on three dimen- 
sional ordering in spin- 1/2 chain compounds where the 
exchange interaction along the chains J is much stronger 
than the inter-chain coupling J 7 . Before doping, the 
Hamiltonian is 
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where J is the site-index along the chains and S are the 
vectors to neighboring chains. The chain lattice spac- 
ing has been set to unity. Randomly removing some of 
the spins breaks the chains up into finite segments with 
open boundary conditions, which are still weakly coupled 
to neighboring chains. This model describes Zn doped 
Sr2Cu03, for example. For the pure system a standard 
method to study the Neel temperature for weakly coupled 
chains is to first determine the staggered susceptibility of 
the one-dimensional chains, %i (T). If we then treat the 



inter-chain couplings in mean field theory [EJ we obtain 
the condition which determines the Neel temperature 



zJ'xi{T N ) = 1, 
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where z is the number of neighboring chains from the sum 
over 6 in Eq. (|l|). Since at low T, Xi(T) diverges as 1/T, 
this predicts T/v oc J', so that Neel order is predicted to 
occur for arbitrarily weak inter-chain coupling. 

In this paper we extend this approach to the doped sys- 
tem by calculating the staggered susceptibility of chains 
with arbitrary length L to find an average value of \i as a 
function of temperature T . This type of mean field treat- 
ment of inter-plane couplings in doped samples was used 
to study the Neel temperature of quasi-two-dimensional 
antiferromagnets in J3|. This method can be carried out 
much more accurately in the quasi-one-dimensional case 
studied here because an analytic expression for the stag- 
gered susceptibility of finite chains can be found, which 
by itself yields rather interesting results, exhibiting very 
different behavior for even and odd length chains. Eq. (||) 
corresponds to approximating the interchain interactions 
as simply providing a staggered field of fixed magnitude, 
acting on a given chain. This approximation results 
from averaging over both quantum fluctuations and im- 
purity locations on neighboring chains. We expect it to 
become more reliable when the average chain length L 
and z (i.e. the lattice dimension) increases. Clearly this 
approximation must break down at large impurity dop- 
ing p = l/L, before the percolation threshold is reached 
(p c sa 69.8% for a three dimensional simple cubic lattice). 
In lower dimensions this approach becomes more ques- 
tionable since the percolation threshold is reached earlier 
and the number of neighbors is lower, while in higher di- 
mensions this method may become exact as z — * oo and 
valid for all doping levels since p c — » 1. 

The staggered susceptibility per unit length of a finite 
chain of length L with open boundary conditions at finite 
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temperature T = 1//3 is 
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In order to find a reliable average we need to determine \i 
for a large range of temperatures and lengths for which we 
use both bosonization techniques and numerical Monte 
Carlo simulations. Note that here we measure the stag- 
gered response to a staggered field, not to be confused 
with the staggered response to a uniform field 0. 

In order to calculate \i as a function of temperature 
we now go to the continuum limit and use the field theory 
treatment which is correct in the asymptotic low T, large 
L limit [Q. The spin operators are then described in 
terms of a boson field, d>, 
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where C is a parameter which will be discussed in more 
detail below. The boson field </> is described by a free 
massless relativistic Hamiltonian |7| up to a marginally 
irrelevant interaction which gives rise to logarithmic cor- 
rections as we will see later. We normalize the operator 
cos-\/27rc6 so that its T = 0, L — oo, equal time correla- 
tion function decays with distance as l/2|r|. In the sum 
of Eq. (||) we can neglect all rapidly oscillating parts so 
that only the second term in Eq. (Q) will be kept in the 
alternating spin-spin correlation function 



G(x, y, t) = C 2 (cos V2tt(/)(x, it) cos \2m(f>(y, 0)) , (5) 
so that Eq. (||) becomes 
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Xi ~ 7 / dT / dx I dyG(x,y,r). (6) 
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To calculate the correlation function we use the mode ex- 
pansion of the boson for a finite chain with open bound- 
ary conditions M 



(j)(x,t) = 4/5 + V2^S z y + <j> > {x,t), 
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contains the ordinary boson modes a n and the "zero 
mode" eigenvalue S z corresponds to the z-component of 
the total spin of the chain, which takes integer values for 
even lengths L and half-integer values for odd L. The 
spin wave velocity is given by v — t:J/2. 

Before considering the case of arbitrary T and L, it is 
interesting to consider the limit T — > with L held fixed. 
In this limit, upon inserting a complete set of states, 



Xi — ► ■j^|(0|S'f lt |0)| 2 will be dominated by the ground- 
state |0), where 5f lt = EjL 1 (-l) J ' +1 S'/. Using Eq. (g) 
we can directly find the local staggered magnetization of 
the lowest energy state in any sector with a given S z 
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For L even S z = in the groundstate so this gives zero 
LT X i -> 0. However, for L odd S z = ±1/2 in the 
groundstates, giving 
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which upon integrating over x gives S z lt oc \fh and there- 
fore xi This divergence is in sharp contrast to 
the even case. Interestingly Eq. ( |Io[ ) indicates a max- 
imum response in the center of the chain which agrees 
with numerical results [fTE] , |l6| and is reminiscent of the 
square-root increase of the staggered response to a uni- 
form field with the distance from the open ends iH. It 
is interesting to note that our finding S z lt oc V ' L for a 
spin chain corresponds to an intermediate result between 
a Neel state with S z lt = ±L/2 and a nearest neighbor 
dimer state with one unpaired spin S z lt — ±1/2. 

We now consider the case of general L and T using the 
field theory approach. The correlation function can then 
be written as 
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Upon using the cumulant theorem for boson modes 
(e A ) = e^ A ^l 2 we can determine the correlation func- 
tion at any finite temperature and length by following 
the analogous calculations in Refs. Q and |l^]. Using 



the shorthand notation u 
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where 6\ is the elliptic theta function of the first kind [[Ll| . 
The parameter 7 = ^7 gives the spacing in the finite 
size energy spectrum in relation to the temperature. The 
contribution B(x) from the zero modes is given by 
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FIG. 1: The scaled staggered susceptibility in Eq. (|1J) 
f(LT) = XiT/C 2 determined from Eqs. (|) and (|l|) for even 
and odd length chains. The upper boundary of the graph 
represents the asymptotic limit as LT —* oo in Eq. @. The 
points correspond to numerical Monte Carlo results for dif- 
ferent combinations of L and T divided by the logarithmic 
correction in Eq. ( |li| ) with a — 23. 



where 8 is the elliptic theta function of the second kind for 
odd chains 8 = 82, while it is the elliptic theta function 
of the third kind for even chains 8 — 03. Remarkably, the 
correlation functions in the continuum limit in Eq. (flT 



therefore retain information about the underlying lattice 
and explicitly depend on the parity of L. This result 
requires the explicit use of the zero modes in the mode 
expansion ||, 
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The difference arises because of the 
different set of eigenvalues of S z : integer and half- integer 
for even and odd length chains, respectively. 

At this point we may rescale all the variables of inte- 
gration in Eq. (0) by L (or alternatively T) to express 
Xi in terms of a universal function of the dimensionless 
variable LT /v. 



Xi = C 2 f(LT/v)/T 



(14) 



In the thermodynamic limit LT/v — > 00, we can use 
the asymptotic behavior of the 9 functions jll| as was 
done in a related calculation in H [l0| for e 1 — ► 1. In 
this case we can combine the two terms of G(x,y,r) in 
Eq. ( |l2| ) into one, giving the known finite temperature 
correlation functions Q. This results in the well-known 
1/T behavior flf] 
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As expected there is no difference between even and odd 
length chains in the limit of LT jv — > 00 in Eq. (|l5|). 

In the opposite limit of zero temperature and finite 
length LT/v — > 0, however, we find a qualitative differ- 
ence for even and odd chains. Again using asymptotic 
limits of 8 functions as e~ 7 — ► 0, we now find 
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Since for odd L the correlation function G(x, y, r) ap- 
proaches a constant as r — > 00 we get a divergence of xi 
in Eq. (||) with 1/T for low temperatures, while for even 
L the integral is proportional to L resulting in 
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L even 



0.913893 C 2 /T, 



L odd 
(17) 

where E is the elliptic integral of the second kind [ pr| , 
which can also be derived from Eq. (|To|) . Note that the 
scaling behavior with 1/T in the two limits LT/v <C 1 for 
odd L in Eq. @ and LT/v > 1 in Eq. (jT||) is the same, 
up to a factor of about 2. This is of crucial significance 
for the behavior of the Neel temperature of the doped 
quasi-one-dimensional system as we shall see. 

So far we have ignored the marginally irrelevant inter- 



action mentioned earlier. Its effect on the staggered sus- 
ceptibility at finite T but L — > 00 is well-known. It corre- 
sponds to replacing the constant C 2 by a slowly varying 
function of T 



C 2 = 20n(aJ/T)/(27r) 3 , 



(18) 



where a is a dimensionless constant |12], |T^|. From fit- 
ting our susceptibility data, we find a ~ 23, which gives 
C 2 w 0.33 ± 0.07 for 0.001 < T/J < 0.2. A finite 
length L together with open boundary conditions leads 
to more complicated logarithmic corrections which may 
in general involve a different exponent near the bound- 
ary [Q. Hence, the general expression of the logarithmic 
correction at finite L is not known, but we expect that 
C 2 m 0.33 ± 0.07 remains approximately correct for the 
relevant length scales studied here. 
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FIG. 2: The Neel temperature Tjv/J' as a function of J' 
and average chains length L from bosonization results (lines) 
compared to Monte Carlo simulations (symbols). 



the curve largely unchanged for different J', and only 
shifts it up by a few percent as J' is lowered. Therefore 
we can make a nearly universal quantitative prediction 
for all doping levels and coupling strengths. 

The Neel temperature is strongly affected by doping 
when the impurity concentration, l/L> J 1 J J and may 
drop by as much as a factor of 5, although it remains 
finite. This is because the scaled average staggered sus- 
ceptibility gCTjfL) of odd chains is finite as L — > so 
that Eq. (09) can always be fulfilled for a positive Tjv. 
If, however, only even chains were allowed in the system, 
no non-zero solution would exist for L < Q.&J/zJ' . As 
mentioned above, we expect this result to break down 
at larger impurity doping as the percolation threshold is 
approached and Neel order disappears. 
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The full behavior of \i as a function of the scaling 
variable LT is shown in Fig. (|l|) compared to numerical 
Monte Carlo data after dividing by the logarithmic factor 
in Eq. (fl8|). There are no adjustable parameters for this 
fitting except for the constant inside the logarithm a ~ 23 
and all numerical points from the Monte Carlo simula- 
tions fall close to this universal line for larger values of 
LT as well (not shown). The errors are less than the size 
of the symbols in the figure so that the deviations are due 
to higher order corrections. For the simulations we chose 
different values of 20 < L < 120 and 0.025 < T/J < 0.2. 
For LT > 4 J the even and odd cases are virtually indis- 
tinguishable, but as LT — ► there is a clear difference in 
the behavior. 

We now want to determine the Neel temperature by 
using the mean field treatment of inter-chain couplings 
in Eq. (|^) , where \i now represents the one-dimensional 
susceptibility averaged over chain lengths with the prob- 



ability distribution, P(L) 



-L/L 



/L, where L is the 



average chain length, corresponding to an impurity con- 
centration of p — 1/L. Because of the scaling form in 
Eq. ( 14 ) it is straightforward to show that the mean field 
condition in Eq. (H) can be written as 



T N = C 2 zJ'g{T N L) 
where g{TjqL) = J dy ye~ y fiTj^Ly) is the 



(19) 



'(T N Ly) is the average of 
the scaling function / in Fig. (nh and Eq. (14) and C 2 is 
given in Eq. ( [l8| ) with T = Tjy. We sec that the solutions 
for the Neel temperature Tjy/ J' in Eq. ( [l9| ) are functions 
of the scaling variable LJ 1 as shown in Fig. (|J) for z = 4 
compared to the results from Monte Carlo simulations. 
The marginal operator leads to weak logarithmic correc- 
tions to this scaling behavior, which leaves the shape of 
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